Thermodynamics of statistically interacting quantum gas in V dimensions 
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We present the exact thermodynamics (isochores, isotherms, isobars, response functions) of a 
statistically interacting quantum gas in T> dimensions. The results in 2? = 1 are those of the ther- 
modynamic Bethe ansatz for the nonlinear Schrodinger model, a gas with repulsive two-body contact 
potential. In all dimensions the ideal boson and fermion gases are recovered in the weak-coupling and 
strong-coupling limits, respectively. For all nonzero couplings ideal fermion gas behavior emerges 
for 2? S> 1 and, in the limit T> —> oo, a phase transition occurs at T > 0. Significant deviations from 
ideal quantum gas behavior are found for intermediate coupling and finite T>. 



I. INTRODUCTION 

The wave of experimental studies that led to the first 
observations of Bose-Einstein condensation (BEG) and 
the development of measurement and confinement tech- 
nologies have renewed strong interest in the statistical 
mechanics of interacting quantum gases P, d, d, 0, H, [1] . 
This line of research can make good use of explicit high- 
accuracy results from any type of analysis that goes be- 
yond low-density /high-temperature expansions and be- 
yond mean-field theory. Of particular interest are results 
for response functions, the very quantities most directly 
amenable to experimental investigations. 

Such results can be produced on a rigorous basis for 
quantum gases with statistical interaction under very 
general circumstances as shown by Wu 0, S| ■ The con- 
cept of statistical interaction introduced by Haldane @ 
has proven to be a very useful methodological device to 
capture the statistical mechanical properties of degrees 
of freedom subject to dynamical interaction. For several 
model systems in dimension V = I the coupling between 
degrees of freedom can be substituted by a generalized 
Pauli principle with no loss of rigor regarding the ther- 
modynamic analysis [ol [lOl. [Til [Ta. fisj . 

Whereas an equivalence between dynamical and sta- 
tistical interaction is not likely to be realized in I? > 1 
(apart from highly contrived scenarios), models of statis- 
tically interacting degrees of freedom can stand on their 
own. Their thermodynamic properties can be analyzed 
exactly in any dimension V, producing a full and consis- 
tent account of fluctuations as will be demonstrated in 
this work [T2| . The exact results emerging from this anal- 
ysis make it possible to connect features of the statistical 
interaction with features of a corresponding dynamical 
interaction. The systematic study of such connections, 
in turn, opens the door to the design of (exactly solv- 
able) models of statistical interaction for the description 
of thermodynamic phenomena associated with specific 
aspects of dynamical interaction. 

In a previous paper [isj we have established a bench- 
mark in that regard by exploring the thermodynamics 
of an ideal quantum gas with fractional statistics in 
V dimensions - a thermodynamic generalization of the 
Calogero-Sutherland (GS) model [M [13, [Hi - taking 



advantage of techniques and results reported in previ- 
ous studiesM,^, 12, 17, IS, M, HH, i22, 23., 25, 2j, 
113, m, [li, |32|, M, Si • In that case the statistical 
interaction was limited to pairs of particles with identical 
momenta. 

Here we relax that constraint and consider a model 
system, again in V dimensions, with a statistical inter- 
action that extends to pairs of particles with arbitrary 
momenta, a system moreover, whose statistical interac- 
tion in 2? = 1 is equivalent to the dynamical interaction 
of a model that is solvable (beyond thermodynamics) via 
Bethe ansatz: the nonlinear Schrodinger (NLS) model 
M, 36, 37,, M, M, |4fl, ii] . 

In Sec. [Ilj we review the concept of statistical interac- 
tion and its use in statistical mechanics. We introduce 
the NLS model and the generalization of its thermody- 
namics to I? > 1. In Sec. mil we describe the method of 
thermodynamic analysis applied to the generalized NLS 
model. In Sec. IIVI we discuss selected thermodynamic 
properties thus calculated. In Sec. |V| we assess the re- 
sults in relation to existing benchmarks for ideal quan- 
tum gases with fractional statistics. 



II. STATISTICAL INTERACTION 

The statistical interaction of any given model system 
is specified by a generalized Pauli principle Q, express- 
ing how the number of states available to one particle is 
affected by the presence of other particles: 



Ad, 



(1) 



The indices i,j refer to particle species and the gij are 
statistical interaction coefficients. For bosons we have 
gij — and for fermions gij ~ 6ij. Integrating Eq. ([T]) 
yields the holding capacity for particles of species i in 
the presence of a specific number of particles from each 
species: 



(2) 



where Ai are statistical capacity constants. The number 
of many-particle states composed of {n^} statistically in- 
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teracting particles is 



(3) 



The three principal specifications of a system of par- 
ticles subject to a statistical interaction are sets of (i) 
energies e.^, (ii) capacity constants Aj, and (iii) interac- 
tion coefficients gij . The grand potential of such a system 
can be expressed in the form ^] 



n = -fesT^AJn 



1+Wi 



Wi 



(4) 



where the Wi are determined by the nonlinear algebraic 
equations, 



kRT 



ln(l + Wi) - ^gji In 



(5) 



The control variables are T (temperature) and /i (chem- 
ical potential). The average numbers of particles, (rii), 
are related to the Wi by the linear equations. 



(6) 



If gij — giSij then all Eqs. ([5]) and (O are decoupled 
and the statistical interaction reduces to a (fractional) 
exclusion condition. 



A. Application to quantum gas 

For a nonrelativistic quantum gas in a box of dimen- 
sionality V and volume V = LP the aforementioned spec- 
ifications are encoded in the energy-momentum relation 
eo(fc) = |kp (in units where /2m = 1) and in a function 
f;(k, k'). The grand potential @ becomes 

n^-ksT(If] (d-kln^-±^^ (7) 
\2ttJ J 

where is the solution of the nonlinear integral equa- 
tion 

=ln(l + ii;k)- / d^fc'g(k',k)lnii^. (8) 



w-k' 



The particle density in k-space, (rik), is the solution of 
the linear integral equation 

{n^)w^+ I dPk'g{]^,\^)M = l. (9) 



The fundamental thermodynamic relations (thermody- 
namic and caloric equations of state) depend on the so- 
lutions of IHl) and (IHl) as follows: 



ksT V27r 
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=(^\ fd-klnl±^^ (10) 



(11) 



(12) 



If the statistical interaction is of the form gdk — k'|) then 
the solutions of ((8]) and ^ only depend on the magnitude 
of the particle momenta. 



B. Nonlinear Schrodinger model 

Consider the boson gas in I? = 1 with repulsive con- 
tact interaction of strength c as described by the NLS 
Hamiltonian 



H 



(13) 



The thermodynamic Bethe ansatz (TBA) solution [33, 
38, 39, 40] of the NLS model expresses the grand poten- 
tial in the form 

17(r, L, ^,) = -ksT 1^°° dk In (l + W A-^) , 

where e(fc) is the solution of the Yang- Yang equation [37| . 
keT 



with kernel 



27r 



K{k~k') 



dk'Kik-k') In (l + e-<'^')/kBT 

(15) 



2c 



c2 + (fc - k'f 



(16) 



The particle density (rifc) is the solution, for given e(fc), 
of the Lieb-Liniger equation [H, [s^l , 



{nk) 



l^(,<k)/kBT 



= 1 + ^ / dk' K{k^k'){nk'). 



(17) 

Bernard and Wu showed that this TBA solution is 
equivalent to the thermodynamics of a statistically inter- 
acting gas in I? = 1 if the following identifications are 
made: 



Wk 



t(k)/kBT 



j{k - k') = d{k - k') - Y^^'^ ^ 



Generalization of NLS model 



(18a) 
(18b) 



The generalized NLS model is a quantum gas in T> 
dimensions with the statistical interaction expressed by 
the kernel 



Wk 



K{k-k') ^ 



7r^/2-ir(P/2) [c2 + (k - k')2]^ 



(19) 
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of the Yang- Yang equation p5)) and Lieb-Liniger equa- 
tion p?)) generalized toV >1 and designed to reproduce 
the exact thermodynamics of the dynamically interacting 
NLS model in I? = 1 . The kernel (fTO)) has the properties 



lim K{k) = 0, 

C— i-OO 

limii:(k) = 27r(5(k), 
y (FkK{k) = 27r. 



(20a) 
(20b) 

(20c) 



This model interpolates between the ideal Fermi-Dirac 
(FD) gas in the strong-coupling limit (c = oo) and the 
ideal Bose-Einstein (BE) gas in the weak-coupling limit 
(c = 0) in all dimensions T> > 1. In the limit I? ^ oo it 
turns into the ideal FD gas for all c > 0. 

For the further analysis of the generalized NLS model 
we reduce Eqs. ([8]) and ([9]) into integral equations for the 
functions 

e(fc) = fcsTlnwk, n{k) = (nk), (21) 
where k = |k|. We also introduce scaled quantities 
k _ . c 



k = 



knT' 



(22a) 



n{k) = n{k). (22b) 



Equations p^)) and ([TT]) thus thermodynamically gener- 
alized to I? > 1 become 

e(k) = P -Inz^ / dfc'i^(fc,fc')ln(l + e-''(^')), (23) 
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FIG. 1: Iterative solutions (a) e{k) + lnz of Eq. ^3^, and (b) 
n(fc) of Eq. (|24p for c = 0.25, z = 3, and various values of D. 
The curves for D — oo are solutions of Eqs. (|29|) and (|30|l . 
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= 1+ / dk'K{k,k')n{k'), (24) 



n{k) 



with fugacity z = e^^^'^^ and reduced kernel 



^.r 7n _ 2c^r(I?) 
K{k,k)- pp/2)j^ 



+ 2c2(P+fc'2) + (P-fc'2)2](I'+l)/2- 



(25) 



Any particular solution e(fc), n{k) at fixed z, c describes 
the system over a range of temperature T, chemical 
potential /x, and coupling constant c. For couplings 
< c < oo the solutions of Eqs. (pS)) . (|24p also depend 
on P. High-precision data for e{k), n{k) can be obtained 
from an iterative solution. The sample of data shown in 
Fig. [1] exhibit the main characteristic features of these 
functions. 

The function e{k) -I- Inz is monotonically increasing 
from a minimum at fc = 0. The dependence on z, c 
(not shown), and T) is strongest at small k and weakens 
rapidly at fc > 1. The presence of a shoulder combined 



with a deep minimum at fc = signals a BEC ordering 
tendency (weakly evident here in the V > 2 data). The 
function n{k) is monotonically decreasing from a maxi- 
mum at fc = 0. Again there is a significant dependence 
on z, c, and V at small fc. The function itself vanishes 
rapidly at fc > 1. 

The dependence on fc, /i, T of the functions e(fc) and 
n(fc) is not reducible, for < c < oo, to a dependence on 
a single variable, (fc^ — nj/ksT. This represents a major 
deviation from a signature property of ideal gases. 



D. Limiting cases 

In the strong-coupling and weak-coupling limits, which 
are fixed points of the scaling relation for the coupling 
constant, we recover the familiar FD and BE results: 



e(fc) — k — Inz, n(k) — \ —e" 

' z 



(26) 
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for c = oo, and 
'1 



e{k) = In i^-e'' - ij , n{k) = l^-e" - ij (27) 

for c = 0. The ^-independence of these functions is an- 
other signature property of ideal gases, a property upheld 
in the presence of fractional statistics [15| . The parabolic 
curve, e(fc) -I- Inz = fc^, may serve as baseline for the so- 
lutions of (|23p. It is exact in the strong-coupling limit 
for all z > 0, but only for z — s- in the weak-coupling 
limit. For < z < 1 the deviations of e{k) -I- Inz from fc^ 
and of fi{k) from zero are suppressed by factors ~ e~'^ 
as k increases. We have noted very similar behavior in 
the numerical results for < c < oo. 

In the limit V ^ oo the reduced kernel ([^5)) acquires 
a much simpler structure. 



hm K 



(k, k') :^S(k' ^ VWT^^ 



(28) 



but the c-dependence is retained and (for c > 0) also the 
statistical coupling of particles with distinct momenta. 
Equations (gS]) and with reduced kernel (UHl) turn 

into the implicit functions 



e{k) = fc2 - In z - In 1 



1 + n ( V fc^ + 



n(fc) l-fe''^*^^ 
The solution of Eq. ^ is 

OO 

and the solution of Eq. ([50)1 is 

oo 

n(fc)=^a;z'+ie-^-'('+i 



(30) 



(31) 



(32) 



1=0 



with the recursion 



ai = e 2 



ai-jC 2 



■J 0-1) 



3 = 1 



-%-(('-j)('-i+i)+2i) 



(33) 



For c — > oo Eqs. (p6)) are recovered, and for c ^ 
Eqs. dSZl). 



III. THERMODYNAMIC ANALYSIS OF 
GENERALIZED NLS MODEL 

Exact results for the thermodynamics of the general- 
ized NLS model in P > 1 and across the range < c < oo 
of coupling strengths can now be calculated from the so- 
lutions e{k) and n{k) of Eqs. and (I24p . respectively. 



A. NLS functions 

The fundamental thermodynamic relations l|10p-(fT2 
are rewritten in the form 



V 



V 

T _ p(V) 



(z,c) 



1 - z 



where 



At — 



/l2 



/2m=l 



2'KmkBT 



(34) 



(35) 



(36) 



(37) 



is the thermal wavelength and the term in psp enclosed 
by square-brackets is relevant only if c = and T) > 2. 
The NLS functions in (l34l)-(l36l) are defined as follows: 



(29) F(^)(z,c)- 



r(x>/2) 



dk k 



ln(l + e-^"W) , (38) 



r(2?/2) Jo 



dkk^-^n{k), (39) 
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Fu\^^^)- + d'k'k-^^'m, (40) 

where the additional dependence of e(fc), n{k) on V, c, z 
is implied. 

In the strong-coupling and weak-coupling limits, the 
NLS functions turn into the familiar FD functions, 



fn{z) 

and BE functions, 



r{n) Jo z-ie^ + 



z>0, (41) 



gn{z) 



1 r°° dxx"-! 



y, 0<z<l, (42) 



r(n) Jq z ^e^ 
respectively: 

Fp(^)(z,(X3) = F^^\z,^) = /p/2+i(z), (43a) 

(43b) 



F^^(z,oo) = /p/2(z). 



and 



Fp(^)(z, 0) = F^^)(z, 0) = .gi,/2+i(^), (44a) 
F^^)(z,0)=52./2(z). (44b) 
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Furthermore, for 2? ^ 1 fermionic behavior results for 
any c > 0: 



F, 



AT 



(45a) 
(45b) 



With increasing I? the factor k^^^ pushes aU significant 
contributions to the integrals ((38l) - (|40)) toward larger and 
larger k, where the deviations of e(k) and n(k) from their 
(c = oo)-values become smaller and smaller [42 1. 

A characteristic ideal-gas property is that the depen- 
dence of the fugacity on the thermodynamic variables 
T, V,Af is expressible as a function of a single variable, 



X = Xtv-^/^, V = V/Af. 



(46) 



In the Maxwell-Boltzniann (MB) gas we have x''^ = z, 
in the FD gas — fv/2{z), and in the BE gas x^ — 
9v/2{z)- A unique functional relation persists in the case 
of fractional statistics [l^ . In the generalized NLS model. 



however, we have a;'^ = Fj^ {z,c) with a separate T- 
dependence contained in c. For ideal quantum gases, 
including those with fractional statistics, there also exists 
a unique functional dependence oipV/MkBT on z. Again 
this no longer holds in the generalized NLS model, where 



we have pV/NkBT = F^'^\z, c)/F^>{z, c) 



CT 



(T = const), 



^JkBTp 




{p — const) 



where x and y are the solutions of 

X {z^CyX')^ 

y^+-^^F(^\z,c,yl 



(51) 



(52) 



(53) 
(54) 



respectively. 

Reference values based on the chemical potential 
present themselves as an alternative in some situations. 
Defining Inz = Ty/T in isochoric processes and Inz = 
Tp/T in isobaric processes, we have 



Ty^ _ PV_ ^ 

Tv Pv 

for V — const, and 



V 
2" 



V 

2^ 



r 



V 



+ 2 



1 



(55) 



(56) 



for p = const. The divergence of these ratios of in the 
limit V oo has some surprising consequence as will be 
discussed in Sec. IIV Bl 



B. Reference values 



IV. RESULTS 



We introduce reference values for the thermodynamic 
variables v, T, p based on the thermal wavelength At and 
the MB equation of state pv — ksT in the presentation 
of our data below: 



47r 



47r 



VT 



knTn 



477 

k^ 



V/2 



Pv 



PT 



,;2/V+l 



{kBT) 



const. 



(47) 



U/2+1 



4.(f 
47r 



(477)^^/2 



(T const.) 

(48) 

{p = const.) 



(49) 

They are especially useful in comparative plots that en- 
compass the full range of c at finite V. 

For the thermodynamic analysis we must adapt the 
NLS functions to the type of process under consideration. 
Each function has a different z-dependence at fixed c, de- 
pending, for example, on whether we consider v — const, 
T = const, or p = const. To this end we introduce three 
kinds of reduced coupling constants for use in isochoric, 
isothermal, and isobaric processes, respectively: 



y/ksTy 



c 

X 



(v = const), (50) 



In Ref. [15| we presented a panoramic view of the ther- 
modynamics of the generalized CS model (ideal quan- 
tum gas with fractional statistics) in T) dimensions. The 
emphasis was on the crossover between boson-like and 
fermion-like features in isochores, isotherms, isobars, re- 
sponse functions, and the speed of sound as caused by as- 
pects of the statistical interaction that reflect long-range 
attraction and short-range repulsion. 

The generalized NLS model considered here exhibits 
some similarities with the generalized CS model regard- 
ing thermodynamic properties, especially their depen- 
dence on the coupling constants of the two models. How- 
ever, there are notable differences, many of which can 
be identified as significant deviations from ideal gas be- 
havior. In our presentation of results we highlight these 
deviations and the role of dimensionality. 



A. Isochores, isobars, and isotherms 

The dependences of p on T at w = const, of u on T at 
p — const, and of p on u at T = const are determined by 
P4p and psp in parametric representations. 



P_ 

Pv 



Fp (Z, CyX^ 



[Fj:^\z,Cyx)] 



T 



1+2. ' T 



F')^\z,Cyx) 



(57) 
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V 

Vr, 



T 



(58) 
(59) 



respectively, with the fugacity z in the role of the param- 
eter. Here x and y are the solutions of ([55)1 and ([5l)) . 
respectively. 

In Fig. [2] we show isochores, isobars, and isotherms for 
various coupling strengths c^^pj- in 2? = 3. The vari- 
ation of the curves between the (weak-coupling) boson 
limit and the (strong-coupling) fermion limit is similar 
to what was observed in the generalized CS model [isj : 
the convergence of all curves toward the MB line at high 
T or large v, and the fanning out at low T or small v. 
Corresponding plots in other 2? show similar trends in 
the two models. 

The shape of the curves for Cy^p^r > in Fig. [5] yield 
some insight into the physical interpretation of the sta- 
tistical interaction. For weak couplings {cy^p^T < 1) the 
curves exhibit boson-like features at high T or large v 
and fermion- like features at low T or small v. These 
observations translate into a long-range attractive part 
and a shorter-range repulsive part of the statistical in- 
teraction. The attractive tail is only present for small 
Cv,p,T, whereas the repulsive core is conspicuous for all 
Cv,p,T > 0. 

Among all the curves only the ones pertaining to the 
boson limit {cy_p^T = 0) have a singularity. This singular- 
ity signals the presence of a phase transition, the onset of 
BEC. In the T> — 3 case shown, the phase transition oc- 
curs at Tc/Ty ~ 0.527, Pc/Pv — 0.271 along the isochore, 
at Tc/Tp ~ 0.889, Vc/vp ~ 0.456 along the isobar, and at 
Vc/vt — 0.383, Pc/pt — 1-341 along the isotherm. 

The bosonic isochore has a singularity at Tc/Ty > 
only in I? > 2. In 2 < I? < 4 it has a discontinuity 
in curvature. In I? > 4 it becomes a discontinuity in 
slope. In the limit T) oo the bosonic isochore itself 
becomes discontinuous. By contrast, the bosonic isobar 
has a singularity at Tc/Tp > in all dimensions V > I, 
but with Vc/vp > only in I? > 2. The bosonic isotherm 
has a horizontal portion at v/vt < Vc/vt in I? > 2 (see 
Ref. [l^ for more details on the bosonic curves.) 

We have already noted that all three NLS functions 
([55)l - (|40p converge toward the corresponding FD func- 
tions as I? — > oo provided we have c > 0. One reflection of 
this fact in the data for isochores, isobars, and isotherms 
is that all curves for c^^p.T > move closer together as 
V increases. They coalesce into the universal curve (iso- 
chore, isobar, or isotherm) representing the ideal FD gas 
in I? = cxD. Only the bosonic curves at T < Tc or v < Vc 
stay apart. 

In Fig. [3] we show two sets of isochores, isobars, and 
isotherms for the generalized NLS model in I? = 1,2,3, 
one set for weak coupling, the other for strong coupling. 
Also shown (dashed) are the corresponding curves per- 
taining to D = oo, which will be derived in Sec. IIVBI 
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FIG. 2: Isochores, isobars, and isotherms in 7? = 3 for Cv,p,T 
0, 0.1, 0.25, 0.5, 1, 3, oo (from bottom up in each frame). 



For the most part, the weak-coupling and strong-coupling 
curves are located on opposite sides of the dashed line in 
each frame. 

Convergence of the data for I? = 1,2,3 toward the line 
representing 2? = oo is only apparent at high T or large 
V and more clearly in the strong-coupling data than in 
the weak-coupling data. This is not surprising in view 
of the observation made earlier in the context of Fig. [2] 
that for weak couplings the (effectively) long-range at- 
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FIG. 3: Isochores, isobars, and isotherms in D = 1,2,3 for 
Cv,p,T = 0.25 and Cv,p.T ~ 3. The dashed hnes represent the 
curves for © = 00. 



B. Phase transition in 17 = 00 

It is well known that no phase transition at T > ex- 
ists for free fermions in I? < 00. No transition is expected 
to exist in the generalized NLS model in finite V except 
in the boson limit. However, a curious transition does 
emerge in the limit V 00, where the generalized NLS 
model with c > effectively turns into an ideal FD gas. 

To determine the thermodynamic equation of state, 
f{p, V, T) — 0, of the generalized NLS model in I? = 00 
we recall (gS]) and rewrite (I571)-(ISni) for I? > 1 with the 
FD functions substituted for the NLS functions. A sin- 
gularity at T > results as a consequence of the fact 
that the two limits I? — > 00, z — s- cx) are not interchange- 
able. The emergence of the singularity is apparent in the 
isochores and isobars in I? ^ 1 as shown in Fig. 21 




FIG. 4: Isochores (main plot) and isobars (inset) of the ideal 
FD gas in D = 1, 2, 3, 10, 40, 300 (top to bottom). The result 
for I? = cxD is represented by the dashed lines. 

The isochore at z < cx) in the limit I? — > 00 yields a 
straight-line segment with unit slope and zero intercept 
in the (T/Ti,, p/p„)-plane over a nonzero interval < 
T < 00: 



P_ •D>1 T fv/2+l{z) P^o c T 

PV Ty fv/2{z) Ty 



(60) 



The reference values (|T7)) become ksTy — An, py — Att/v 
in the limit I? — > 00. The same isochore in the limit V — > 
00, z —^ 00 with 'D/2 = rlnz, r > yields a horizontal 
line segment over a nonzero interval < ^ < J^c: 



tractive part and short-range repulsive part of the statis- 
tical interaction are responsible for opposing trends and a 
crossover between them. However, convergence becomes 
manifest in higher V (not shown) as the NLS functions 
gradually turn into FD functions first for strong couplings 
and then also for weak couplings. 



p P»i fv/2+i{z) 



1-1-2/2?' 



(61) 



[fv/2{z)] - 2 h^' 



where we have used the leading term in the asymptotic 
expansion of the FD function |4^. The value of Tc is 
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determined by the intersection point of the two Hne seg- 
ments. The equation of state thus reads 44 1 



pv 



knTr 



An 

e 



(63) 



This same universal relation can also be inferred from 
Eq. ([58|) for the isobar or from Eq. ((59|) for the isotherm 
by performing the appropriate limits. All isotherms are 
hyperbolas, including the transition line at T = Tc- All 
the isochores and isobars consist of two straight-line seg- 
ments with the singularity at T = T^, as already shown 
in Fig. m 

This somewhat unusual phase transition from a fully 
intact Fermi sea at T < Tc to an ideal MB gas at T > Tc 
results from the conspiracy of two opposing effects, one 
suppressing thermal excitations at low T and the other 
enhancing them at high T. Both effects grow stronger in 
higher dimensions. 

We know from that as V increases the reference 
temperature T„ becomes smaller and smaller compared 
to the Fermi temperature T^ in isochoric processes (con- 
sidered here for specificity). This suppresses any rise in 
pressure at sufficiently small but nonzero T/Ty more and 
more strongly. In the limit T> —^ oo, as T^/T^ — *■ 0, the 
pressure will remain constant over a non-vanishing in- 
terval of T/Ty at the value p/pv — e^^ exerted by the 
perfect Fermi sea. 

We also know (e.g. from analogies to microcanonical 
ensembles) that as T> increases the energy density of one- 
particle states is progressively thinned out inside the sur- 
face of the Fermi hypersphere except close to the surface. 
The consequence is that a smaller and smaller amount of 
thermal energy is needed to knock out the vast majority 
of particles from the Fermi sea. Moreover, the density of 
vacancies near the Fermi edge becomes so large that the 
constraint on occupancy imposed by the Pauli principle 
is negligible. 

In I? ^ 1, therefore, if T/Ty is raised gradually, no sig- 
nificant thermal excitations take place initially because 
Ty/Ty ^ 1. The isochore stays flat. Once T has reached 
a certain threshold the Fermi sea is emptied quickly be- 
cause of its shallowness and the abundance of vacancies 
close by. The system thus crosses over from a near per- 
fectly degenerate Fermi sea to a nearly ideal MB gas on 
a very short interval of T/Ty as documented in Fig. 3) In 
T> = oo this crossover has sharpened into a phase transi- 
tion. There is no latent heat involved in that transition 
and there is no sudden increase in pressure [45| . 

Note that on the alternative temperature scale Ty the 
emergent crossover between near perfect Fermi sea and 
almost ideal MB gas is pushed to lower and lower values 
of T as I? increases, ultimately to T/Ty = for I? = oo. 
The resultant isochore is then that of the ideal MB gas 
all the way down. 

It is interesting to recall the phase diagram of the ideal 
BE gas in 2? = oo for comparison. The thermodynamic 
equation of state inferred from the scaled isochores, iso- 
bars, or isotherms as derived, for example, in Ref. (l5| 



has the form 



pv 



:T, T >Tc 
T <Tc 



knTc 



47r. 



(64) 



As in the FD case there are two phases separated by a 
transition line at constant T. The high-T phase is again 
an ideal MB gas. The low-T phase is a pure BEG. The 
transition is of first order and occurs at a higher temper- 
ature than in the FD case. Whereas the FD transition 
disappears in I? < oo, the BE transition persists down to 
T> > 2, but is of second-order in P < oo and occurs along 
a line in (p, w, T)-space that is no longer an isotherm. 



C. Response functions 

The three major response functions for a gas 
of spinless particles are the isochoric heat capacity, 
Cy = M^^{dU/dT)y, the isobaric expansivity, ap = 
v^^{dv / dT)p, and the isothermal compressibility, kt = 
—v~^{dv/dp)T- For the generalized NLS model we must 
evaluate the expressions 



Cy 



TpUp — 



p2 

4 



Tp 
T 



V 

" 2" 

V 

2" " 



Ptkt 



F, 



F^\z,CyX) 



FPiz,Cpy)-§;F;^-'iz,Cpy) 



(65) 
V 



Fj:^\z,Cpy)^F^''\z,Cpy) 2 



p(I'), 



V ^F^\z,ct) 

VT f^F^-^\z,CTy 



(66) 



(67) 



versus the independent variables T/Ty, T /Tp, v/vt, re- 
spectively, from ((57)1 - ((59)) . 

In Fig. [5] we show the dependence of each response 
function on the coupling strength in I? = 3. The 
variation of the curves between the BE and FD limits 
shows some resemblance to that observed in an ideal 
gas with fractional statistics (generalized CS model) p^ . 
All three response functions depend only weakly on the 
statistical interaction at high temperature or low den- 
sity. The dominant trends there reflect MB behavior, 
Cy — (3/2)fcsr, ap = l/T, kt = v/ksT. Distinct 
boson-like and fermion-like features and crossovers be- 
tween them emerge at low temperatures and high densi- 
ties. The exact analytic behavior of the response func- 
tions in any V for the FD and BE limits was described 
in Ref. [15]. 

The heat capacity in I? = 3 for strong coupling is 
dominated by fermion-like features at all T, exhibiting a 
monotonic descent from the MB asymptote as T is low- 
ered and a linear approach to zero. For weak coupling 
the initial increase, the smooth maximum followed by a 
steep descent is a boson-like feature. The ultimate lin- 
ear approach to zero signals the crossover to fermion-like 
behavior. 
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FIG. 5: Isochoric heat capacity, isobaric expansivity, and 
isothermal compressibility in D = 3 for Cv,p,T = 0, 0.1, 0.25, 
0.5, 1, 3, cxD. The heat-capacity curves for c„ = 3, oo are 
unresolved. The dashed lines mark the onset of BEG. 



The expansivity in I? = 3 depends only weakly on T for 
strong coupling and approaches zero linearly as T ^ 0, 
which is a characteristic fermion-like behavior. For weak 
coupling the pronounced rise in expansivity is a boson- 
like feature. However, the repulsive core of the statistical 
interaction for c > 0, no matter how weak, prevents the 
expansivity from diverging and forces the fermion-like 
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FIG. 6: Isochoric heat capacity in D = 1,2 for Cv ~ 0, 0.1, 
0.25, 0.5, 1, 3, oo. The dashed lines reflect the high-T asymp- 
totic values Cv/ks = V/2. 



behavior at low T. 

Stiff resistance to compression, perhaps the most out- 
standing fermion-like feature, makes itself manifest with 
growing strength in the strong-coupling compressibility 
curves in I? = 3 as the density is increased. In the weak- 
coupling curves, on the other hand, we observe trends 
reminiscent of bosonic behavior at moderate densities. 
While the BE curve diverges at v = Vc, the repulsive 
core of the statistical interaction prevents the transition 
from taking place if c > 0. The compressibility curve 
bends into a smooth maximum or a mere shoulder down 
to fermionic stiffness. 

Similarities to the response functions of an ideal quan- 
tum gas with fractional statistics are also manifest in 
other dimensions V. However, there are two notable ex- 
ceptions. In Fig. [Hlwe show the heat capacity for I? = 1 , 2 
in the same format as the data for 2? = 3 in Fig. \5\ 

InV — 2 the isochoric heat capacity of an ideal quan- 
tum gas is well-known not to depend on the exclusion 
statistics [11,124, 25, 26, H [Ml, [13 . That is no longer the 
case in the presence of a statistical interaction such as re- 
alized in the generalized NLS model. Only the two curves 
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representing the weak-coupling and strong-coupling lim- 
its coincide. The curves at intermediate coupling are 
subject to shifting trends caused by the long-range at- 
tractive and short-range repulsive parts of the statistical 
interaction. 

In I? = 1 the heat capacity curves are monotonically 
increasing functions near the BE limit and functions with 
one smooth maximum near the FD limit. Upon variation 
of the exclusion statistical parameter between the two 
limits in the quantum ideal gas the appearance of the 
smooth maximum coincides with a switch in sign of the 
leading correction to the high-T asymptote [15|. Upon 
variation of the statistical coupling strength between the 
same limits in the generalized NLS model, the smooth 
maximum at low T appears before the approach to the 
asymptote switches side from below the asymptote to 
above it. In consequence there is a range of coupling 
strengths where the heat capacity has a smooth maxi- 
mum followed by a smooth minimum. 

In Fig. [7] the focus is on the P-dependence of the three 
response functions at c^,p.T — 0.25 (weak coupling) and 
Cv,p,T — 3 (strong coupling). The goal is to gain further 
insight into how gases with boson-like, fermion-like, and 
crossover features in 2? = 1, 2, 3 evolve into one and the 
same FD system as 2? — s- oo with an emergent singu- 
larity at T/Ty = T/Tp — e^^ in isochoric and isobaric 
processes. 

For the scaled heat capacity we show weak-coupling 
results in the main plot and strong-coupling results in the 
inset, the dependence on coupling strength of the results 
in I? = 1, 2, 3 is conspicuous but becomes imperceptibly 
small for V > 20. The universal FD result for I? = oo 
(dashed line) is a simple step function. 



lim - — 



= e{T-Tc), Tc/Ty = e- 



(68) 



The decrease in initial slope with increasing V is clearly 
visible in the curves for V = 1,2,3 but overall conver- 
gence toward the step function is slow. The scaled heat 
capacity of the ideal BE gas in I? = cxd also consists of 
two terms, one being a step function as in (|68p but with 
Tc/Ty = 1 and the other being a i5-function representing 
the latent heat [isj . 

For the isobaric expansivity the weak-coupling and 
strong-coupling curves are shown in the same plot. There 
is very little variation between V ^ 1 and P = 3 for 
the strong-coupling case {cp = 3). Somewhat larger and 
more systematic variation occurs in the weak-coupling 
case (cp = 0.25). The universal FD result for V = oo 
(shown dashed) is 
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FIG. 7: Isochoric heat capacity (scaled), isobaric expansivity, 
and isothermal compressibility in D = 1, 2, 3 for reduced cou- 
pling strengths c^,ji,t = 0.25 and Cv,p,T ~ 3. The dashed lines 
represents the FD curves for D = oo. To avoid cluttering we 
have omitted one expansivity curve. 



Tr 



Tc/Tp 



(69) 



Convergence is slow, but evident in the curves for V = 
40, 300. Comparing the resuh ([Ml) for the ideal FD gas 
with that of the ideal BE gas, both in I? = oo, we find 
that the latter also has the form (IM)) at T > Tc but with 



Tc/Tp = 1. The expansivity of the FD gas is zero at 
T < Tc, in the BE gas it is undefined [15]. 

For the isothermal compressibility the universal FD 
line for I? = oo (dashed line), 



PtKt 



VT 



(70) 
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is indistinguishable from the MB result. The weak- 
coupling and strong-coupling curves are located on op- 
posite sides of that line. Convergence is apparent in 
the curves for "D = 1,2,3 in the strong-coupling case 
(ct = 3) but not in the weak-coupling case {ct ~ 0.25). 
The isothermal compressibility of ideal BE gas is also de- 
scribed by the result (|70p but only for v/vt > Vq/vt — 1- 
At v/vt < ^ the bosonic result is infinite (15j |. 



D. Speed of sound 

The speed of sound as inferred from c = (pKs)~^/^, 
where p = m/v is the mass density and ns the adiabatic 
compressibility, can be brought into the form [46j 



mc 



{v/vt) 



{T/T,)\v/vT){T,a,)^ 

{PTKT){Cv/kB) 



(71) 



which only involves dimensionless quantities previously 
determined in terms of the NLS functions, 



knT 



I 2 + F< 
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Here c must be replaced by c^, c^x^ or c^y depending on 
whether we are considering an isothermal, isochoric, or 
isobaric process, respectively. In Fig. [5] we present data 
for the T-dependence of the speed of sound of a weak- 
coupling system under isobaric and isochoric conditions. 

It is well-known that in ideal gases the curves for 
m(? jkBT differ from those of the isobars or isochores 
only by a multiplicative factor (1 -t- 2/1?). We have seen 
that this relation still holds in the presence of fractional 
statistics [l^. The data shown here for the generalized 
NLS model demonstrate that no such relation holds any 
longer in the presence of a statistical interaction that is 
not reducible to a simple exclusion principle. 

The deviations appear to be strongest in I? = 1. 
All deviations are expected to fade away in the limt 
I? — > 00 when ideal gas behavior is restored as explained 
in Sec. IIVBI Particularly noteworthy is the observation 
that the T-dependence of the speed of sound in I? = 1 
at constant (average) pressure undergoes a minimum as 
highlighted in the inset. No such minimum exists in the 
isobar. 

We attribute this effect to the crossover between boson- 
like features at high T and fermion-like features at low T . 
The general trend, realized in ideal gases is that the speed 
of sound decreases monotonically upon cooling. Super- 
imposed on this is another trend that signals softness 
when boson-like features are predominant and stiffness 
when fermion-like features are predominant. 



FIG. 8: (a) Main plot: Speed of sound (squared and scaled), 
mc? /ksTp, at constant (average) pressure versus T /Tp in 
T> — 1,2,3 (solid lines) and D — 00 (dotted line). Also shown 
are scaled isobars, {l + 2/'D){v/vp) versus T/Tp (dashed lines). 
Inset: Zoomed extract of the T> = 1 results near the mini- 
mum of the speed of sound data, (b) Corresponding data for 
mc^/ksTy and {l + 2/T>){p/p^) versus T/T^. All data pertain 
to a weak-coupling situation {cv,p = 0.25). 



For stronger coupling (e.g. c„.p = 3) the deviations 
of the speed-of-sound data from the scaled isobar or iso- 
chorc are of a similar kind and size. We have detected 
no minimum at T > in these data. If such a minimum 
exists at all it must occur at very low T, out of reach of 
our numerical analysis. 



V. CONCLUSION 

The exact thermodynamic analysis of the generalized 
NLS model, a quantum gas in V dimensions with a sta- 
tistical two-body interaction, has yielded significant de- 
viations from characteristic ideal quantum gas behavior 
in several respects. 

For given coupling strength < c < cx) (i) the aver- 
age level occupancy {n{k)) is no longer a unique function 
of (fc^ — fij/ksT and independent of "D; (ii) the quan- 
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titles pX^/ksT, MX^/V, and {U\l^/V)/{kBTV/2) are 
no longer unique functions of the fugaclty z; (111) the 
two quantities pX^/ksT and {U\^ /V) / {ksTV /2) are 
no longer Identical. 

Among the consequences are (1) that the T-dependence 
of the Internal energy Is no longer of the same shape as 
the Isochore; (11) that the quantity pV^/A/'fesT Is no longer 
a function of z alone In given T); (111) that there Is no 
longer any simple relation between the speed of sound 
and the Isochore or Isobar. 

In any finite T) the statistical Interaction of the gener- 
alized NLS model smoothly Interpolates between an Ideal 
BE gas In the weak-coupling limit (c = 0) and an Ideal 
FD gas In the strong-coupling limit (c = oo). In the 
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